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Abstract 

We show that in see-saw models with bimaximal lepton mixing at the GUT scale and 
with zero CP phases, the solar mixing angle #12 generically evolves towards sizably 
smaller values due to Renormalization Group effects, whereas the evolution of #13 
and 023 is comparatively small. The currently favored LMA solution of the solar 
neutrino problem can thus be obtained in a natural way from bimaximal mixing 
at the GUT scale. We present numerical examples for the evolution of the leptonic 
mixing angles in the Standard Model and the MSSM, in which the current best-fit 
values of the LMA mixing angles are produced. These include a case where the 
mass eigenstates corresponding to the solar mass squared difference have opposite 
CP parity. 
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1 Introduction 



Recent experimental evidence strongly favors the LMA solution of the solar 
neutrino problem with a large but non-maximal value of the solar mixing angle 
0i2 [1-4]. An overview of the current allowed regions for the mixing angles and 
the mass squared differences is given in table 1. 





Best-fit value 


Range (for 6 


G [0°,45°]) 


C.L. 


012 [°] 


32.9 


26.1 - 


43.3 


99% (3a) 


023 [°] 


45.0 


33.2 - 


45.0 


99% (3a) 


013 [°] 




0.0- 


9.2 


90% (2a) 


Am 2 sol [eV 2 ] 


5 • 10~ 5 


2.3 • 10~ 5 - 


3.7- 10~ 4 


99% (3a) 


lAmLJ [eV 2 ] 


2.5 • 10~ 3 


1.2 • 10~ 3 - 


- 5 • 10~ 3 


99% (3a) 



Table 1 

Experimental data for the neutrino mixing angles and mass squared differences. 
For the solar angle 6*12 and the solar mass squared difference, the LMA solution has 
been assumed. The results stem from the analysis of the recent SNO data [3], the 
Super-Kamiokande atmospheric data [5] and the CHOOZ experiment [6]. 

A big problem for model builders is to explain the deviation of 12 from max- 
imal mixing, while keeping 823 maximal and #13 small at the same time. The 
Renormalization Group (RG) evolution is a possible candidate for accomplish- 
ing this. Therefore, it is interesting to investigate the evolution of the mixing 
angles from the GUT scale to the electroweak (EW) or SUSY-breaking scale. 
A number of studies with three neutrinos considered the possibility of increas- 
ing a small mixing angle via RG evolution [7-10]. Others focused on the case 
of nearly degenerate neutrinos [11-16], on the existence of fixed points [17], or 
on the effect of non-zero Majorana phases on the stability of the RG evolution 
[18]. 

We consider the see-saw scenario, i.e. the Standard Model (SM) or MSSM 
extended by 3 heavy neutrinos that are singlets under the SM gauge groups 
and have large explicit (Majorana) masses with a non-degenerate spectrum. 
Due to this non-degeneracy, one has to use several effective theories, with the 
singlets partly integrated out, when studying the evolution of the effective 
mass matrix of the light neutrinos [19,20]. Below the lowest mass threshold, 
the neutrino mass matrix is given by the effective dimension 5 neutrino mass 
operator in the SM or MSSM. The relevant RGE's were derived in [20-25]. 

In this paper, we assume bimaximal mixing at the GUT scale with vanishing 
CP phases and positive mass eigenvalues. We calculate the RG running nu- 
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merically in order to obtain the mixing angles at low energy and to compare 
them with the experimentally favored values. We include the regions above 
and between the see-saw scales in our study, which have not been considered 
in most of the previous works. We find that the solar mixing angle changes 
considerably while the evolution of the other angles is comparatively small, 
so that values compatible with the LMA solution can be obtained. We present 
analytic approximations that help to understand this behavior and show that 
it is rather generic. 



2 Bimaximal Mixing at the GUT Scale 



At the GUT scale, we assume bimaximal mixing in the lepton sector. We 
restrict ourselves to the case of positive mass eigenvalues and real parameters, 
so that there is no CP violation. In the basis where the charged lepton Yukawa 
matrix is diagonal, up to phase conventions the general parametrization of the 
effective Majorana mass matrix of the light neutrinos is then 



m 



bimax 



V{\ , 0, f ) • diagK, m 2 , m 3 ) • V T {% 0, f ) 



a — o c —c 



cab 
—cba 



(1) 



where 



V(9i2, #13, #23) 



C12C13 S12C13 S13 

-C23S12 - S23S13C12 C23C12 - S23S13S12 S23C13 
S23S12 - C23S13C12 -S23C12 - C23S13S12 C23C13 



(2) 



with s 



sin 6ij and = cos 6^- is the (orthogonal) CKM matrix in standard 



parametrization, and 



a = - (mi + m 2 + 2 m 3 ) , 



- (-mi - m 2 + 2m 3 ) , 
m 2 — mi 
2V2 ' 



Inverting equations (3) yields the mass eigenvalues 

mi = a — b — V2c , 
m 2 = a — b + V2 c , 
m 3 = a + b . 



(3a) 
(3b) 
(3c) 



(4a) 
(4b) 
(4c) 
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From equation (3a) we see that a > 0. Equations (4) imply that the solar mass 
squared difference Am^ ol = m\ — m\ is related to c, while the atmospheric one, 
Am^ tm = mj — mf , is controlled by b. Thus, a > \b\ > \c\. For b > we obtain 
a normal mass hierarchy, while for b < the mass hierarchy is inverted, as 
illustrated in figure 1. For positive c, ni\ < rri2, otherwise mi > m%. Hence, 
Amg 0l is positive only if c is. If a ^> |c|, the spectrum is called degenerate. 
We use the convention that the mass label 2 is attached in such a way that 
< #i2 < 45°. This can always be accomplished by the replacement c <-> — c. 

m 3 m 2 /i 

m 2 /i 

m l/2 ffl3 



(a) Normal hierarchy (b) Inverted hierarchy 

Fig. 1. Possible mass hierarchies for the light neutrinos. We use the convention that 
mi and m,2 are chosen in such a way that < 612 < 45°. The LMA solution then 
requires ni2 > mi. 



In our see-saw scenario, the effective mass matrix of the light neutrinos is 

mf™ = ^Y?M- 1 Y v (5) 

at the high-energy scale, with (4>) = w 174 GeV. Obviously, the singlet 
Yukawa and mass matrices Y v and M cannot be determined uniquely from 
this relation, i.e. there is a set of \Y V , M} configurations that yield bimaximal 
mixing. After choosing an initial condition for Y v , M (and thus the see-saw 
scales) is fixed by the see-saw formula (5) if Y v is invertible. 



3 Solving the RGE's 



To study the RG running of the leptonic mixing angles and neutrino masses, 
all parameters of the theory have to be evolved from the GUT scale to the 
EW or SUSY-breaking scale, respectively. Since the heavy singlets have to be 
integrated out at their mass thresholds, which are non-degenerate in general, 
a series of effective theories has to be used. The derivation of the RGE's and 
the method for dealing with these effective theories are given in [20]. Starting 
at the GUT scale, the strategy is to successively solve the systems of coupled 
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differential equations of the form 



d («) (n) / r (n) \ 

167rV^ = ^AUi}J (6) 

(») («) r (n) (n) (n) -1 

for all the parameters X i} Xj G \K, Y U ,M, . . . / of the theory in the energy 
ranges corresponding to the effective theories denoted by (n). At each see- 
saw scale, tree-level matching is performed. Due to the complicated structure 
of the set of differential equations, the exact solution can only be obtained 
numerically. However, to understand certain features of the RG evolution, an 
analytic approximation at the GUT scale will be derived in section 5. 



4 Examples for the Running of the Mixing Angles 



Figures 2 and 3 show typical numerical examples for the running of the mixing 
angles from the GUT scale to the EW or SUSY-breaking scale. They contain 
an important effect that appears for most choices of the initial parameters: 
The solar angle 812 changes drastically, while the changes in #13 and 823 are 
comparatively small. This agrees remarkably well with the experimentally fa- 
vored scenario. 




Fig. 2. RG evolution of the mixing angles from the GUT scale to the SUSY-breaking 
scale (taken to be ~ 1 TeV) in the MSSM extended by heavy singlets for a normal 
mass hierarchy and Y u = X diag (l, s, e 2 ) with tan/? = 5, e = 0.525, a = 0.0675 eV 
and X = 1. In this example, the lightest neutrino has a mass of 0.025 eV. The 
kinks in the plots correspond to the mass thresholds at the see-saw scales. The 
grey-shaded regions mark the various effective theories. 
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Fig. 3. Example for the RG evolution of the mixing angles in the SM extended by 
heavy singlets from the GUT scale to the EW scale for a normal mass hierarchy and 
Y v = X diag (l,e,e 2 ) with e = 0.65, a = 0.0655 eV and X = 1. In this example, 
the lightest neutrino has a mass of 0.024 eV. 

5 Analytic Approximation for the Running of the Mixing Angles 
at the GUT Scale 



In order to understand the effect found numerically in the previous section, 
we now derive an analytic approximation for the RG evolution of the mixing 
angles at the GUT scale. It is only affected by the part of the RGE that is not 
proportional to the unit matrix, which is given by 



167T /J, 



d/i 



-m,. 



a 



YX] m v + C e m v \Y}Y ( 



C v 



yW„ 



m v + C v m v 



Y}Y„ 



+ terms with trivial flavour structure 



(7) 



with C P 



2> 



G v = \ in the SM and C e = C u = 1 in the MSSM. Analogously 
to equation (1), m v is parametrized by 

m v (t) = v(e 12 (t),9 13 (t),e 23 (t)) ■ m diag (t) • v T (e 12 (t),e 13 (t),e 23 (t)) , (8) 



where /i is the renormalization scale, t := In and mdi ag := diag(m 1; m 2 , m 3 ) 



In general, the real Y v can be written as 

Y* = V(£ 12 , ^13, 62) • diag(j/i, y 2 , y 3 ) ■ V T ((f) 12 , 13 , 32 ) . 



(9) 



However, the effective mixing matrix m u is invariant under the transformations 
Y v — > V T Y V and M — > V T MV, which correspond to a change of basis for the 
heavy sterile neutrinos. Thus, V(£\ 2 , £13, £32) in equation (9) can be absorbed 
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into M, leading to the simpler parametrization 



Y v (y u y 2 , ya, 0i2, 0i3, 032 ) = diag(yi, y 2 , y 3 ) • V T (0i 2 , 0i 3 , ^32) 



(10) 



Furthermore, we use the approximation that the effect of the charged lepton 
Yukawa matrices Y e can be neglected compared to that of the neutrino Yukawa 
matrix. Note that in the MSSM a large tan/5 can yield a relatively large Y e , 
which can also have sizable effects that are neglected in this approximation. 

We now differentiate equation (8) w.r.t. t and insert the RGE (7). For the 
evolution of the mixing angles at the GUT scale with bimaximal mixing as 
initial condition, we thus obtain both in the SM and in the MSSM the ratios 



012 
013 



2 \J2 (mi + 777,2) (^3 — m±) (^3 — m 2 ) Fi 



m gut K m 2 ~ mi , 



8 (m! ~~ m l m 2) F 2 + 4 V2 (777 2 — 777i) 772.3 F 3 



± — = — — — - — r for hierarchical neutrino masses 1 



'2^/2 m 2 - mi F 2 
1 Am' F! 



11a) 



for degenerate neutrino masses 



012 
023 



I 2v/2 Am s 2 ol F 2 

2 (mi + m 2 ) (m 3 - mi) (m 3 - m 2 ) Fi 



m gut (m 2 -mi) 8 (m 2 - m x ) m 3 F 2 + 4 V2 (mj - mi m 2 ) F 3 

(lib) 



± — — - — i- for hierarchical neutrino masses 1 

2m 2 -m 1 F 3 

1 A< m Fi 



I 2 Am s 2 ol ^3 



for degenerate neutrino masses 



with 



Fi 



(cos 



,(2 13 ) - 3) sin(2 23 ) - 6 cos 2 (0i 3 ) 

(2 023)) 



(yl-y 2 2 ) {cos(2 

-4 cos(2 23 ) sin(2 0i 2 ) sin(0i 3 )} 

+ (yl + y 2 2 -2 yf) [ cos(2 13 ) (sin(2 23 ) - 3) + (1 + sin(2 23 )) " , 

(12a) 



2 ivl-yt) cos(0i 3 ) sin(2 0i 2 ) (sin(0 23 ) - cos(0 23 )) 
- (y 2 + y 2 2 - 2 y 3 2 + (y\ - yj) cos(2 i2 )) x 
x sin(2 0i 3 ) (cos(0 23 ) + sin(0 23 )) , 



(12b) 



Note that this approximation is also valid for a relatively weak hierarchy, where 
m 3 is a few times larger or smaller than mi, m 2 . 
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^3 = (vl-yf) [cos(2 12 ) (cos(2 13 ) -3) cos(2 23 ) 
+ 4 sin(2 0i 2 ) sin(0i 3 ) sin(2 23 ) 
+ 2 (yl + yl -2yl) cos 2 (0i 3 ) cos(2 23 ) . 



(12c) 



This result can also be obtained from the formulae derived in [26]. The con- 
stants Fi, F 2 and F 3 clearly depend on the choice of Y^(Mqut)- However, 
unless the parameters {y±, y 2 , y 3 , 0i 2 , 0i 3 , 32 } are fine-tuned, we expect the 
ratios Fi/F 2 and Fi/F 3 to be of the order one. Consequently, the RG change 
of 9 12 is larger than that of the other angles if the mass-dependent factors 
in equations (11) are large. This is always the case for degenerate neutrino 
masses, since Am^ tm ^> Aro s 2 ol . As (mi— m 2 ) is related to the small solar mass 
squared difference, it is also true for non-degenerate mass schemes, unless mi 
is very small, in which case the ratio approaches 1. This corresponds to a nor- 
mal mass hierarchy and a strongly hierarchical mass scheme. Finally, it can 
be shown that the running of 9yi is always enhanced compared to that of 6>i 3 
and # 23 for inverted schemes. Hence, we conclude that this is a generic effect. 

6 Parameter Space Regions Compatible with the LMA Solution 

6.1 Parameters at the GUT Scale 

The considerable change of the solar mixing angle found in the previous sec- 
tions raises the question whether the parameter region of the LMA solution 
might be reached by RG evolution, if one starts with bimaximal mixing at high 
energy. We will investigate this possibility by further numerical calculations 
in the following. To reduce the parameter space for the numerical analysis, we 
choose a specific neutrino Yukawa coupling Y v at the GUT scale. We assume 
that it is diagonal and of the form 



Y v and M are now determined by the parameters {e,X, a, b, c}. Moreover, 
we fix the GUT scale values of b and c by the requirement that the solar 
and atmospheric mass squared differences obtained at the EW scale after the 
RG evolution be compatible with the allowed experimental regions. Thus, 
we are left with the free parameters X, e and a. The parameter e controls 
the hierarchy of the entries in Y v and thus the degeneracy of the see-saw 
scales, while a determines the mass of the lightest neutrino. The dependence 
of physical quantities on e and a is shown in figure 4. The effect of changing 
the scale X of the neutrino Yukawa coupling will be discussed in section 6.3. 
As mentioned above, we work in the basis where the Yukawa matrix of the 
charged leptons is diagonal. 




(13) 
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0.04 



0.15 

a[eV] 



0.25 



0.10 



0.55 



0.99 



e 



(a) Parameter a 



(b) Parameter e 



Fig. 4. Plot 4(a) shows the mass of the lightest neutrino (at low energy) as a 
function of a for the SM and the MSSM with normal mass hierarchy, X = 1 and 
e £ [0.1,0.99] (grey region). Plot 4(b) shows the degeneracy of the see-saw scales, 
parametrized by \n{M^/Mi) (at the GUT scale), as a function of e for the same 
cases with a G [0.04 eV, 0.25 eV] (grey region). 

6.2 Allowed Parameter Space Regions 

The parameter space regions in which the RG evolution produces low-energy 
values compatible with the LMA solution are shown in figure 5 for the SM 
and the MSSM (tan/3 = 5) with a normal mass hierarchy. We find that for 
the form of Y u under consideration, hierarchical and degenerate neutrino mass 
schemes as well as degenerate and non-degenerate see-saw scales are possible. 
For inverted neutrino mass spectra, allowed parameter space regions exist as 
well. 

We would like to stress that the shape of the allowed parameter space regions 
strongly depends on the choice of the initial value of Y u at the GUT scale. 
One also has to ensure that the sign of Am^ is positive, as the LMA solution 
requires this if the convention is used that the solar mixing angle is smaller 
than 45°. With bimaximal mixing at the GUT scale, the sign of Am^, is not 
defined by the initial conditions. Using the analytic approximation of section 
5, the sign just below the GUT scale can be calculated. We find Am^ > 
for F\ < and vice versa. However, in order to predict the sign of Am^ 
at low energy, the numerical RG evolution has to be used. This excludes 
some of the possible choices for the neutrino Yukawa coupling Y v at the GUT 
scale. For example, among the possibilities with diagonal Y v it excludes Y u = 
diag(e 2 ,e, 1). 
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0.04 0.15 0.25 0.04 0.15 

a [eV] a [eV] 



0.25 



(a) SM 



(b) MSSM 



Fig. 5. Parameter space regions compatible with the LMA solution of the solar 
neutrino problem for the example Y v = diag(l, e, e 2 ). The initial condition at the 
GUT scale M GU T = 10 16 GeV is bimaximal mixing, and the comparison with the 
experimental data is performed at the EW scale or at 1 TeV for the SM and the 
MSSM, respectively. The white regions of the plots are excluded by the data (LMA) 
at 3<r. For this example, we consider the case of a normal neutrino mass hierarchy 
and X = 1 for the scale factor of the neutrino Yukawa couplings. 

6. 3 Dependence on the Scale X of the Neutrino Yukawa Coupling 

For small values of X, the contribution from Y v to the evolution of the mixing 
angles above the largest see-saw scale is suppressed by a factor of X 2 . Never- 
theless, the evolution to the LMA solution is still possible, as can be seen from 
the example in figure 6. Here the large change of 9 12 also seems to be generic 
but takes place between the see-saw scales, which shows the importance of 
carefully studying the RG behavior in these intermediate regions [20]. Note 
that in this case the analytic approximation of section 5 cannot be applied, 
since it is only valid at the GUT scale. 

6.4 Effect of Neutrino CP Parities 

An example for the running of the mixing angles to the LMA solution with a 
negative CP parity for the state with mass m 2 is shown in figure 7. For this 
we have chosen a different diagonal structure for Y u , 



at the GUT scale. Here, the evolution to the LMA solution is possible due to 
running between the see-saw scales. A more detailed study of the effect of CP 
phases will be given in a forthcoming paper [27]. 
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Fig. 6. RG evolution in the SM for X = 0.01, e = 0.3, a = 0.0535 eV and a normal 
mass hierarchy. The running from bimaximal mixing to the LMA solution now takes 
place between the see-saw scales. In this example, the lightest neutrino has a mass 
of 0.017 eV. 
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Fig. 7. RG evolution in the SM with a negative CP parity for m,2, X = 0.5, 
e = 3.5 • 10~ 3 and a normal mass hierarchy. The running from bimaximal mixing 
to the LMA solution takes place between the see-saw scales. In this example, we 
consider a strongly hierarchical mass spectrum. The lightest neutrino has a mass 
of 0.004 eV. Note that the cases fl 23 > 45°, Am| 3 > and 9 23 := 90° - fl 23 < 45°, 
Am-23 := — Am-23 are indistinguishable in neutrino oscillations. 

The large RG effects in this case seem surprising at first sight, since previous 
studies, e.g. [18,26], found that opposite CP parities for mi and m 2 prevent 
a sizable change of the solar mixing angle by RG evolution. However, these 
works did not consider the energy region between the see-saw scales, where 
the largest change occurs in our example. This fact explains the apparent 
discrepancy. 
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6. 5 Low Scale Values of #13 and 623 

The mixing angles 9 13 and 9 23 are affected by the RG evolution as well, i.e. 
they do not stay at their initial values #13 = 0° and 623 = 45°. However, lower 
bounds on their changes cannot be given unless a specific model is chosen. As 
one can see from the previous examples, the changes can be tiny. For instance, 
the evolution of figure 2 gives A9 13 = 0.02°, which corresponds to sin 2 (2# 13 ) = 
5 • KT 7 , and A9 23 = 0.28°. On the other hand, other choices of Y„ at the GUT 
scale produce A#i3 and AQ23 that come close to the experimental bounds. 
This can make it possible to discriminate between models with different initial 
values Y^Mqut). 



7 Summary and Conclusions 

We have shown that in see-saw scenarios the experimentally favored neutrino 
mass parameters with the LMA solution of the solar neutrino problem can be 
obtained in a rather generic way from bimaximal mixing at the GUT scale by 
Renormalization Group running. We have concentrated on the case of vanish- 
ing CP phases, which implies positive mass eigenvalues. In an example where 
the mass eigenstates corresponding to the solar mass squared difference have 
opposite CP parity, we have demonstrated that an evolution towards the LMA 
solution is possible in this case as well. The general case of arbitrary CP phases 
is beyond the scope of this letter and will be studied elsewhere [27]. The mix- 
ing angles evolved down to the electroweak scale show a strong dependence 
on the mass scale of the lightest neutrino, on the degeneracy of the see-saw 
scales, and on the form of the neutrino Yukawa coupling. A generic feature of 
the Renormalization Group evolution is that the solar mixing angle 612 evolves 
towards sizably smaller values, whereas the change of #13 and #23 is compar- 
atively small. In the SM and MSSM, we find extensive regions in parameter 
space which are compatible with the LMA solution for normal and inverted 
neutrino mass hierarchies and for large and small absolute scales of the neu- 
trino Yukawa couplings. Thus, RG running may provide a natural explanation 
for the observed deviation of the LMA mixing angles from bimaximality. 
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